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ABSTRACT 

Context. It is renown that for solar-mass stars planet formation begins in a circumstellar disc. The study of transitional 
objects exhibiting clear signs of evolution in their discs, such as the growth of dust particles and beginning of disc dispersal, is 
fundamental to understanding the processes governing dust-grain coagulation and the onset of planet formation. 
T-H Aims. We attempt to characterise the physical properties of T Chamaeleontis, a transitional T Tauri star showing UXOri-type 
variability, and of its associated disc, and probe the possible effects of disc-clearing processes. 
' Methods. Different spectral diagnostics were examined, based on a rich collection of optical high- and low-resolution spectra. 
The cross-correlation technique was used to determine radial and projected rotational velocities, shape changes of photospheric 
lines were analysed via bisector-method applied to the cross-correlation profile, and the equivalent widths of both the Li I A6708 A 
r^ . photospheric absorption and the most prominent emission lines (e.g., Hq, H/3 and [Ol]6300A) were measured. Variability in 
O I the main emission features was inspected by means of line-profile correlation matrices. Available optical and near-infrared 
JL photometry combined with infrared data from public catalogues was used to construct the spectral energy distribution (SED) 
t| and infer basic stellar and disc properties. 

Results. Remarkable variability on timescale of days in the main emission lines. Ha changing from pure emission to nearly 
photospheric absorption, is correlated with variations in visual extinction of over three magnitudes, while the photospheric 
I I absorption spectrum shows no major changes. The strength of emission in Ha and H/3 is highly variable and well correlated 
with that of the [Ol] lines. The structure of the Ha line-profile also varies on a daily time-span, while the absence of continuum 
veiling suggests very low or no mass accretion. Variations of up to nearly lOkms"^ in the radial velocity of the star are measured 
^ ^ on analogous timescales, but with no apparent periodicity. SED modelling confirms the existence of a gap in the disc. 

Conclusions. Variable circumstellar extinction is inferred to be responsible for the conspicuous variations observed in the stellar 
I continuum flux and for concomitant changes in the emission features by contrast effect. 

Clumpy structures, incorporating large dust grains and orbiting the star within a few tenths of AU, obscure episodically the star 
and, eventually, part of the inner circumstellar zone, while the bulk of the hydrogen-line emitting-zone and outer low-density 
wind region traced by the [Ol] remain unaffected. In agreement with this scenario, the detected radial velocity changes are also 
explainable in terms of clumpy material transiting and partially obscuring the star. 



Key words. Stars: variables: general - Stars: pre-main sequence - Stars: late-type - Stars: individual: T Cha - Stars: circumstellar 
matter - Planetary systems: protoplanetary discs 
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^ 1. Introduction 

T Tauri stars (TTS) are low-mass (M<2Mq) pre-main 

sequence (PMS) stars (Appenzeller & Mundt 1989). They 

Send offprint requests to: schisano i i i • ^ ^ u i • i 

■^■^ are commonly classined into two sub-groups, the classical 

based on observations obtained at the European Southern mmn / mmn \ i • i i i i ; ■ n ; i • i 

Observatory at La SiUa, Chile in program 63.L0045(A); 65.L TTS (cTTSs), which are surrounded by an optically thick 

0089(A); 66.C-0616(A); 67.C-0155(A); 67.C-0155(B); 68.C- ^^^^ ^^^^^ t^^^ ^^"^^^ material, and the weak TTS 

0292(A); 69.C-0207(A); 70.C-0163(A); 073.C-0355(A); 074.A- (wTTSs), presumably representing the final stages of ac - 



9018(A); 075.C-0399(A-F) cretion and disc-clearing processes (Bertout et al. 20071. 



Tables 1, 2 and 4, and Figures 21-32 are only available in Hence, the cTTS-wTTS dichotomy is ascribed to difFer- 
electronic form via http://www.edpsciences.org ent physical processes associated with the evolution of 

Correspondence to: schisano@oacn.inaf.it these young solar- type stars and their circumstellar en- 
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vironmcnt. The equivalent width of the Ha emission is 
commonly used as an empirical criterion to distinguish be- 
tween cTTS and wTTS ( .White fc Basri|[2003J . However, 
due to possible variability, no clean cut can be defined be- 
tween the two subgroups based on the Ha emission alone. 

Besides the traditional cTTS versus wTTS scheme, 
Herbst et al. (19941 introduced an additional class, the 



early-type TTauri stars (eTTSs), from simple analysis of 
the photometric variability. ETTSs are earlier than KO 
and have large-amplitude irregular variations (up to ~3 
mag in V) with no sign of continuum veiling. Within 
this class, those authors also included some Herbig Ae/Be 
stars (HAEBEs) of the UXOri type ( [Herbst et ar]|1999[ ). 
Studying this kind of variability along with stellar mass. 



Bertout ( 2000J concluded that T Tauri stars can display a 
similar phenomenology. Well-studied cases such as RY Lup 
(|Gahm et al. |1989[ ) and AATau ( [Bouvier et al.|2003[ ) sup- 
port this hypothesis. 

How can all these observations be interpreted coher- 
ently in a scenario of star and circumstellar environment 



evolution is still questioned. Grinin ( 1988, ) modelled this 
kind of variability with dusty clumps in Keplerian or- 
bits that temporarily obscure the central star. The rapid 
temporal evolution of such events indicates that the dust 
clouds usually appear in the disc close to the dust sub- 
limation radius, probably in the inner rim region fNattaJ 



et al. 2001 puUemond et al. 20031. This model imphes 



a high inclination angle, while imaging and near-infrared 
interferometry of some of those stars are inconsistent with 
it (see the review by Millan-Gabet et ar]|2007 and refer- 



ences therein) . Another scenario was proposed by Vinkovic 
& Jurkic (2007), in which a dusty outflow creates a halo 



in the inner-disc region, where clumps of dust occasion- 
ally intercept the line of sight and mask the stellar photo- 
sphere. 

T Tauri discs evolve on a timescale of a few million 
years and the circumstellar raw material from which plan- 
ets are constructed undergoes different processes. In the 
inner part of the disc (<10 AU), the gas is accreted onto 
the star or expelled by bipolar jets, while the outer part of 



the disc may be dissipated by photoevaporation (Clarke 
et al.||2001| [Alexander & Armitage 2007). On the other 



hand, dust grains settle in the midplane, grow in size, and 
form planetesimals ( Sicilia-Aguilar et al.|2007 Dullemond 



|fc Dominik 2005] ). As the disc becomes progressively opt: 
cally thin, solid particles migrate inward or outward, de- 
pending on their size, due to the combined action of dif- 
ferent perturbing forces (photophoresis, gravity, radiation 
pressure, and gas drag), until they reach a stability dis- 
tance from the star ( Hermann et al.|2007 ). It remains un- 
clear how planetesimals are formed in circumstellar discs 
and whether the hypothetical dusty clumps in highly vari- 
able TTS are related or not to the evolution of circum- 
stellar dust until the formation of planetesimals. Yet, de- 
tailed studies of the observed properties of objects with in- 
termediate characteristics, in-between those of cTTS and 
wTTS, are expected to provide key constraints on star and 
planet formation scenarios. 



Observations in the infrared (IR) with the Spitzer 
satellite have provided a wealth of information for inter- 
pretating the spectral energy distribution (SED) of young 
stars (see Lada et ar||2006[ Cieza et ar]|2007[ Evans et al 



2008 



and references therein) , allowing the development of 
new classification schemes and scenarios of disc evolution 
(e.g., |Cieza||2008||Merm et al.||2008] . In particular, young 
stellar objects with no IR excess shortward of 10 /im, but 
with significant excess emission at wavelengths longer than 
^24 /im, are interpreted as being in an evolutionary phase 
of disc clearing and are hence classified as transitional ob- 
jects ( [Furlan et al.||2006[ [Najita et al.||2007[ ). Evidence for 
gaps in the disc of transitional objects have also been re- 
ported ( [Brown et al.|2007[|Espaillat et al.|2008| . In one of 
these, LkCa 15, the most likely mechanism for gap-opening 



in the disc seems to be planet formation (Espaillat et al 



2008). Disc clearing may also be due to the presence of 



a companion, as in the case of CoKu Tau/4 (Ireland & 
Kraus[p008| 



Here we present a synoptic study of the spectroscopic 
variability of T Chamaeleontis (hereafter T Cha) , a tran- 



sitional TTauri star (Brown et al. 2007) displaying a 
UXOri- like behaviour. TCha was observed in the course 
of a project aiming to reveal and monitor young spectro- 



scopic binary systems ([Esposito et al 
ar[ [2007| 



2006 Guenther et 



The scheme of the paper is the following. In Sect.|2[ 
we summarise the main observed properties of T Cha as 
known from previous studies. In Sect.|3[we present the ob- 
servations and data reduction, while in Sect.|4|we describe 
the radial velocity and vsinz determinations. In Sects. |5[ 
and[6| we report the analysis performed on the spectrum of 
the object, focusing on the non-photospheric contribution 
and, in particular, on the variability of the most promi- 
nent emission features (e.g.. Ha, H/3, and the [Oi] 6300 A 
lines). In Sect. [s] we analyse the spectral energy distribu- 
tion of TCha, while in Sect.[9| we discuss the results and 
present our interpretation of the object. 

2. Observed properties of T Cha 

T Cha shows strong photometric variability (up to 3 mag- 
nitudes in the V-band) mainly characterised by erratic 
changes ( Hoffmeister|1958| Mauder fc Sosna|1975 jCovino 



et al.[[l992 l, and by a UXOri-like behaviour (i.e., the star 
becomes redder as it fades), with a tight correlation be 
tween brightness and colour ( Alcala et al.|1993 Covino et 



al. 



1996). A perio dicity of 3'*. 2 was also pointed out by 



Mauder & Sosna (19751. Spectroscopic variability in the 



most prominent emission lines was reported by [Gregorio- 
([19921) and Alcala et al.[ ([1993[) 



Hetem et al. 



The PMS nature of T Cha was established by [Alcala[ 



et al. ( 1993 1 on the basis of the lithium criterion combined 



with the Ha emission strength and the presence of strong 



IR excess emission, while the spectral type G8V (Alcala 



et al. 1993) is earlier than that of a typical TTauri star. 



From the strength of the Ha emission in former spectra. 



T Cha had been initially regarded as a wTTS ( Alcala et 
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consistent with the strong variability of the star ( 


Alcala et 


al.||1995 


Geers et al. 


2006 


Gras-Velazques & Ray 


2005 





Moreover, besides the IR excess indicating the presence 
of circumstellar material, the object occasionally exhibits 
forbidden neutral oxygen lines that are not seen in wTTS 



spectra. Brown et al. ( 2007 1 modelled the SED and inter- 
preted it in terms of a transitional disc. 

T Cha is located near the edge of the small cloud 
G300.2-16.9, also known as the Blue Cloud, in the direc- 
tion of the Chamaeleon dark cloud complex (Nehme et al. 
2008). In the Hipparcos Catalogue, a distance of 66j!i[2 PC 
is reported, although the error in the parallax is probably 
far higher because of the strong variability of the star (re- 
ported to vary between 10.4 and 13.4 mag in the _Hp-band 
during the period of observation) . More reliable estimates 
of the distance were obtained from proper motion stud- 



ies by [Frink et al.| (|1998| and by Terranegra et al.| ( |1999) , 
suggesting that T Cha is part of an association of stars at 
nearly 100 pc. Here we adopt the latter value. 

3. Observations and data reduction 

3.1. High-resolution spectroscopy 

Data were acquired at the ESO-La Silla Observatory us- 
ing the echelle spectrograph FEROS (Fiber-fed Extended 
Range Optical Spectrograph), first installed at the 1.5m 
and, since October 2002, then mounted at the 2.2m tele- 
scope. The high resolving power (R~48000) and the wide 
useful spectral range (3600- -9200 A) make FEROS par- 
ticularly suitable for radial-velocity (RV) monitoring and 
spectral-line variability studies. FEROS was operated in 
the Object + Sky configuration, in which a sky spectrum 
is acquired simultaneously with the object through an ad- 
jacent fiber. 

From 1999 to 2006, we obtained a total of 50 FEROS 
spectra with quite different signal-to-noise ratios (S/N), 
due not only to the change of telescope or different ob- 
serving conditions but also to the strong variability of the 
star. The spectra were acquired within the framework of 
a project on young spectroscopic binary systems. For that 
reason, the time coverage of our data is rather uneven, 
with some spectra acquired on daily temporal base and 
others separated by several months. 

The data reduction was performed using the specific 
FEROS Data Reduction Software (DRS) implemented in 
the ESO-MIDA^ environment. The reduction steps were 
the following: bias subtraction and bad-column masking; 
definition of the echelle orders on fiat-field frames; sub- 
traction of the background diffuse light; order extraction; 
order by order flat-fielding; determination of wavelength- 
dispersion solution by means of ThAr calibration-lamp 
exposures; rebinning to a linear wavelength-scale (AA = 
0.03 A) with barycentric correction; and merging of the 
echelle orders. All the spectra were then normalized to 
the continuum. 



More details about the data reduction procedure and 
technical specifications of the instrument can be found at 
the FEROS Web sitc0 We emphasise that the high sta- 
bility of the instrument allows an internal accuracy in the 
wavelength calibration of approximately 200 ms~^ (aver- 
age residual r.m.s. ~ 3.5 x 10~^ A). A log of the FEROS 
observations is provided in Table [l] 

Apart from the FEROS data, we also report measure- 
ments from a previous set of six high-resolution CASPEC 
spectra obtained simultaneously with optical broad-band 
photometry on two consecutive nights, 1994 January 31 
and February 1 (Covino et al. 1996). Details about the 
instrumental set-up and data reduction can be found in 
Covino et al. (1997). 

3.2. Low-resolution spectroscopy 

In addition to the high-resolution spectroscopy, we in- 
clude unpublished low-resolution (R~2000) spectra gath- 
ered during various observing runs conducted between 
1993 and 1995 using the Boiler & Chivens spectrograph 
at the ESO 1.5m telescope at La Silla, Chile. The reduc- 
tion of these data was performed as described in [Alcala et| 
al. (1995). We recall that these spectra were calibrated in 



relative fluxes using a standard star. 

4. Cross-correlation function analysis 

4.1. Radial velocity determinations 

We used the cross-correlation technique to determine 
the radial velocity of TCha from the FEROS spec- 
tra. As a template, we chose a spectrum of the G8.5V 
star HD152391 (RV=44.8kms~^ and usini = 3kms-i) 



(Nordstrom et al. 20041, observed with the same instru- 



mental set-up. 

We followed the prescription by Esposito et al. (20061, 
i.e., a Gaussian was fltted to the peak of the cross- 
correlation function (CCF), computed in six distinct spec- 
tral ranges. In some cases, the Gaussian flt failed, and in 
such cases we calculated the moments of the CCF profile 
and took the first-order moment as an RV estimator. In 
the end, we adopted as RV value and associated error, 
the average of the RV determined in the different ranges 
and their r.m.s. dispersion, respectively. All RV measure- 
ments are reported in Table [l] The RV values vary with 
time ranging between 6 and 30kms^^. A period search 
was performed on these data, but no clear indication of 
periodicity was found. 

4.2. CCF bisector analysis 

The Gaussian fit may sometimes appear inadequate due 
to asymmetries in the CCF peak. Asymmetries in stellar 
absorption lines can arise for several reasons and may be 
conveniently represented by line bisectors (Gray 2005). In 



Munich Image Data Analysis System 



http: / /www.ls.eso.org/lasilla/sciops /2p2 /E2p2M /FEROS 



4 Schisano et al.: Variability of TCha 

Table 1. Journal of FEROS observations and radial velocity measurements. Only in electronic form. 



Date 




HJD 


RVg 


O'RV 


Date 


HJD 


RV, 


O'RV 








-2400000. 


(km/s) 


(km/s) 






-240000. 


(km/s) 


(km/s) 


27 


March 


1999 


51265.537 


14.27 


0.27 


22 


May 2000 


51687.648 


22.97 


1.77 


27 


March 


1999 


51265.576 


16.89 


0.68 


23 


May 2000 


51688.640 


9.82 


0.31 


28 


March 


1999 


51265.662 


15.09 


1.63 


04 


January 2001 


51913.860 


13.86 


2.13 


28 


March 


1999 


51266.573 


18.21 


0.73 


09 


January 2001 


51918.775 


15.65 


2.15 


29 


March 


1999 


51266.636 


21.39 


0.68 


14 


January 2001 


51923.824 


14.89 


0.56 


29 


March 


1999 


51266.672 


21.30 


0.90 


19 


April 2001 


52019.578 


16.00 


1.81 


30 


March 


1999 


51267.692 


14.09 


0.25 


26 


April 2001 


52026.586 


14.48 


1.12 


30 


March 


1999 


51267.801 


17.88 


0.47 


26 


April 2001 


52031.573 


14.43 


0.28 


30 


March 


1999 


51267.874 


18.01 


1.38 


14 


February 2002 


52319.745 


12.63 


0.36 


31 


March 


1999 


51268.624 


11.32 


0.87 


17 


February 2002 


52322.742 


27.48 


2.56 


31 


March 


1999 


51268.727 


16.21 


0.86 


08 


April 2002 


52372.606 


14.70 


- 


31 


March 


1999 


51268.850 


6.43 


2.01 


21 


April 2002 


52385.529 


16.44 


0.87 


31 


March 


1999 


51268.900 


12.20 


1.44 


01 


May 2002 


52395.617 


4.33 


1.22 


01 


April 1999 


51269.657 


12.05 


1.15 


02 


May 2002 


52396.614 


9.93 


- 


01 


April 1999 


51269.771 


13.09 


1.00 


12 


March 2003 


52710.601 


21.73 


0.30 


01 


April 1999 


51269.857 


15.80 


1.56 


12 


March 2003 


52717.620 


14.13 


1.51 


16 


May 1999 


51315.622 


18.17 


0.76 


25 


March 2003 


52723.588 


15.96 


2.33 


17 


May 1999 


51316.612 


6.43 


0.82 


02 


April 2004 


53097.541 


15.57 


0.87 


18 


May 1999 


51317.627 


16.67 


1.77 


09 


April 2004 


53104.652 


11.28 


1.18 


19 


May 1999 


51318.591 


13.50 


1.72 


26 


April 2004 


53122.685 


15.93 


0.48 


21 


May 1999 


51320.591 


17.78 


0.52 


10 


May 2004 


53137.672 


22.23 


0.66 


22 


May 1999 


51321.580 


12.72 


2.92 


10 


May 2005 


53501.652 


29.06 


0.87 


09 


May 2000 


51674.670 


18.52 


0.46 


02 


June 2005 


53524.597 


6.63 


1.35 


19 


May 2000 


51684.617 


9.90 


0.52 


30 


June 2005 


53552.546 


14.31 


1.09 


21 


May 2000 


51686.632 


21.19 


2.59 


19 


August 2005 


53602.488 


29.51 


1.95 



some cases, variations in the bisector shape may arise due 
to photospheric spots crossing the stellar disc as the star 
rotates, whereas uncorrelated variations in both bisector 
position and shape may indicate an unresolved spectro- 
scopic system ( [Santos et aL]|2002| |Dall et al.||2006| ). 

From the correlation between RV changes and bisector 
orientation, Queloz et al. ( |2001 j, in the case of HD 166435, 



stated that the periodicity in their RV data was not due to 
an unseen secondary object orbiting the star, but to stel- 
lar activity. We adopted the same tool to verify whether 
RV shifts measured by a Gaussian fit are purely due to 
changes in the CCF peak shape, or reflect true changes 
in the radial velocity of the star. An example of the vari- 
ability of the CCF peak (and corresponding bisector) of 
TCha on daily timescales is shown in Figs. [l]andj2l For 
each of the six spectral ranges considered in Sect. |4~1| we 
divided the corresponding CCF bisector into two inter- 
vals as inlQueloz et al. (20011: an upper part, where the 



strongest changes occur, and a lower one, where the bisec- 
tor position appears more stable (see in Fig. [2]). The RV 
difference between the top and bottom bisector's mean ve- 
locities yields information about the bisector orientation. 
This quantity, called "bisector velocity span" (Toner & 
Gray||1988| [Queloz et al.|2001| [Pall et al.[[2006[ ) , is a good 



measure of the changes in bisector orientation, and hence 
of the asymmetry of the peak. The variations in the cen- 
tral part of the CCF, which reflect the higher stability of 
the wings compared to the peak can be easily discerned. 
The comparison between the RV values obtained from a 



Gaussian-fit and by a bisector analysis is shown in Fig. [3| 
The left panel shows how the RV values from a Gaussian- 
fit are well correlated with the bisector velocity sparj^ 
The Gaussian function reproduces well only one of the two 
sides of an asymmetric CCF peak, leading to a systematic 
shift in the measured RV. On the other hand, the bisector 
values, at the bottom part of the profile, versus the bisec- 
tor span show a weaker correlation than the Gaussian-fit 
values (Fig. |3] right panel). 

The RV determinations derived from the bisector 
method, although of lower amplitude than the values de- 
rived from the Gaussian-fit method, confirm that RV is 
indeed variable, considering the average uncertainty of 
1.7kms~^. We used Fourier analysis, in the formulation 
of the periodogram given by Scargle]( pJ82| , to search for 
possible periodicities in the radial velocity data, although 
no clear values were identified. This might in part be due 
to the extremely uneven temporal sampling. We had only 
one relatively long run of six consecutive nights, but when 
analysing this run alone, no hint of periodicity was found 
on a timescale of a few days. 

Although the presence of blended spectroscopic com- 
ponents cannot be excluded, the lack of a clear periodicity 
prevents us from drawing firm conclusions about the pos- 



sible binarity or multiple-nature of T Cha. As in Queloz 



et al. (2001 1, the correlations between the RV values from 
the Gaussian-fit method and the "velocity span" indicate 



all RV measurements are reported in Tabled available elec- 
tronically 
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Fig. 1. Changes in the CCF in the period 19-23 May 
2000, in two different wavelength ranges. The different 
colours refer to the following dates: 19 May (black), 21 
May (red) 22 May, (green), 23 May (blue). The corre- 
sponding bisectors are also traced. Velocity is relative 
to the template star HD 152391. 



Fig. 2. CCF peaks between 16 and 22 May 1999 (upper 
panel). Velocity is relative to the template star. The 
different colours refer to the following dates: 16 May 
(black), 17 May (red), 18 May (green), 19 May (blue), 
21 May (yellow), 22 May (pink). The lower panel shows, 
on expanded scale, the corresponding CCF bisectors. 



as most plausible explanation the presence of inhomo- 
geneities of variable extent moving across the photospheric 
disc, either intrinsic to the star (e.g., cool spots), or ex- 
ternal to it (e.g., clustered orbiting material). Hence, in 
the following, T Cha will be treated as a single star. For 
the purposes of photospheric spectrum subtraction (see 
Sect. 5.2), we used the bottom bisector's mean RV value. 



4.3. Projected rotational velocity 

We evaluated the projected rotational velocity, wsini, 
of T Cha by measuring the full width at half maximum 
(FWHM) of the CCF peak for the spectra with higher 
signal-to-noise ratio. We first established a relation be- 
tween the FWHM and v sin i, as described in |Covino et al.| 
( 1997 1, using the template spectrum of HD 152391 broad- 



ened artificially by various vsini values. The vsini value 
determined in this way was 37±2kms~^. The same mea- 
surement was repeated on an average spectrum obtained 
by combining a total of 34 spectra with S/N> 50 around 
6000 A. The V sin i determined in the average spectrum 



is consistent with the previous value. However, when ex- 
tending the measurements to all spectra, we identified a 
few with narrower spectral lines, yielding a v sin i of nearly 
30 kms~^. Hence, we checked the v sin i determinations by 
directly fitting each spectrum with an artificially broad- 
ened template, which confirmed that some spectra had 
smaller line widths. 

In Fig. H we report the v sin i values obtained from 
CCF Gaussian-fit FWHM calibration versus the width at 
half maximum of the CCF peak measured directly. In this 
case, it is obvious that the width of the lines does not 
provide a reliable measure of the rotational velocity. A 
variable line-width is indicative of either blended spec- 
troscopic components or the presence of inhomogeneities 
moving across the stellar disc that alter the shape of pho- 
tospheric lines. However, no correlation is found between 
the RV (either from a Gaussian-fit or bisector method) 
and the vsini measurements. 

From a periodogram analysis based on data points 
from the spectra of higher signal-to-noise ratio, we found 
a peak at about 0.25cyclesday~^, or a period close to 
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RV Gaussian (km/s) RV bottom bisector (km/s) 

Fig. 3. The radial velocity measurements from Gaussian-fit method versus the "bisector velocity span" (i.e. the 
difference between the upper and bottom bisector mean velocities) are plotted in the left panel. The right panel 
shows, instead, the mean radial velocity of the CCF bisector bottom versus the "bisector velocity span" . 



35 



50 60 
Width CCF (km/s) 



Fig. 4. Projected rotational velocity from CCF Gaussian- 
fit calibration versus FWHM of the CCF peak. Only the 
data of the highest S/N are considered here. 

4 days. This period does not differ dramatically from that 



found by Mauder & Sosna ( 1975 1 and might reflect a rota- 



tional modulation induced by photospheric spots or other 
types of inhomogeneities transiting over the stellar disc. 

5. The spectrum of T Cha 

5.1. Photospheric spectrum and lithium abundance 

The spectrum of T Cha is that of a G8 V star with strong 
absorption in the Lii resonance line at A 6708 A and a few 



emission lines typical of TTS. The spectral type does not 
change from one spectrum to the other, as also demon- 
strated in Sect. [7) Figure [S] shows the good match be- 
tween the photospheric spectrum of T Cha and that of 
the template. The equivalent width (EW) of the lithium 
line was measured in each spectrum of our FEROS data- 
set. The measurements were found to be internally consis- 
tent, within the errors, yielding a m ean value of EW^i = 
360±30 mA. Using the calibration of Pavlenko & Magazzii 



(1996) and adopting the temperature of 5000 K, we de- 



rived a lithium abundance compatible with the cosmic 
value (log N^i = 3.2), with an uncertainty of 0.2 dex. It 
is important to recall that this determination is strictly 
dependent on the effective temperature. At a temperature 



of 5500 K, as adopted by Kenyon & Hartmann ( 1995 1 for 



a CSV star, the abundance would increase to 3.6 dex, im- 
plying some lithium enrichment with respect to the inter- 
stellar value. 

5.2. Photospheric spectrum subtraction 

One of the main goals of this work was to analyse the 
variability in the emission lines and hence probe the cir- 
cumstcUar environment of TCha. For this purpose, it 
is useful to remove the photospheric contribution. We 
subtracted the the photospheric-line spectrum from the 
FEROS spectra of T Cha using the same template spec- 
trum of HD 152391 used in the cross-correlation analysis 
presented in Sect. |4] 

All spectra were rebinned to the rest-wavelength frame 
in order to account for their different RV shifts, derived 
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Table 2. Radial velocity determinations, both from gaussian fit and bisector method, and equivalent widths of Ha 
and [Oi] 6300A measured on FEROS spectra. Only in electronic form. 



HJD RVg ARVa 


RVb ARVb 




0"Ha 


Wh/3 


0"H/3 


Weaoo 


O"6300 


-2450000 (km/s) (km/s) 


(km/s) (km/s) 


(A) 


(A) 


(A) 


(A) 


(A) 


(A) 



51265.537 


14.27 


0.27 


14.64 


1.07 


-6.26 


1.06 


-1.37 


0.78 


-0.81 


0.05 


51265.576 


16.89 


0.68 


16.54 


1.14 


-6.56 


0.50 


-0.34 


0.29 


-1.00 


0.14 


51265.662 


15.09 


1.63 


17.21 


0.61 


-7.50 


0.30 


-0.56 


0.24 


-1.01 


0.11 


51266.573 


18.21 


0.73 


17.81 


0.54 


-10.64 


0.42 


-2.00 


0.29 


-1.34 


0.01 


51266.636 


21.39 


0.68 


20.57 


1.81 


-12.46 


0.37 


-2.02 


0.38 


-1.48 


0.15 


51266.672 


21.30 


0.90 


20.70 


1.98 


-12.36 


0.80 


-2.89 


0.32 


-1.56 


0.11 


51267.692 


14.09 


0.25 


14.35 


2.16 


-2.89 


0.39 


-0.01 


0.21 


-0.58 


0.12 


51267.801 


17.88 


0.47 


17.58 


0.48 


-2.72 


0.26 


0.03 


0.23 


-0.43 


0.06 


51267.874 


18.01 


1.38 


17.81 


0.82 


-2.55 


0.29 


0.03 


0.22 


-0.34 


0.03 


51268.624 


11.32 


0.87 


15.40 


1.88 


-2.83 


0.36 


0.32 


0.23 


-0.40 


0.03 


51268.727 


16.21 


0.86 


17.10 


1.16 


-5.32 


0.35 


-0.41 


0.20 


-0.39 


0.05 


51268.850 


6.43 


2.01 


18.39 


1.29 


-5.05 


0.41 


0.07 


0.25 


-0.35 


0.04 


51268.900 


12.20 


1.44 


18.75 


2.07 


-3.96 


0.35 


0.10 


0.36 


-0.39 


0.03 


51269.657 


12.05 


1.15 


12.76 


1.46 


-6.18 


0.47 


-0.31 


0.28 


-0.78 


0.05 


51269.771 


13.09 


1.00 


13.42 


0.91 


-5.19 


0.63 


-0.76 


0.40 


-0.77 


0.11 


51269.857 


15.80 


1.56 


16.56 


0.29 


-4.23 


0.55 


-0.95 


0.44 


-0.80 


0.17 


51315.622 


18.17 


0.76 


17.65 


0.83 


-22.68 


0.39 


-2.88 


0.51 


-1.88 


0.25 


51316.612 


6.43 


0.82 


10.54 


3.80 


-23.93 


0.95 


-5.89 


0.70 


-2.06 


0.22 


51317.627 


16.67 


1.77 


20.90 


1.78 


-20.20 


1.03 


-4.30 


0.39 


-0.88 


0.06 


51318.591''^ 


13.50 


1.72 


13.42 


0.49 


-20.65 


0.28 


-3.03 


0.32 


-1.22 


0.07 


51320.591 


17.78 


0.52 


12.61 


2.59 


-1.42 


0.30 


1.08 


0.17 


-0.28 


0.03 


51321.580 


12.72 


2.92 


9.64 


2.72 


-7.56 


0.55 


-0.07 


0.29 


-0.64 


0.04 


51674.670 


18.52 


0.46 


17.39 


1.38 


-0.75 


0.38 


0.87 


0.67 


-0.35 


0.03 


51684.617 


9.90 


0.52 


10.45 


1.69 


-2.39 


0.34 


0.43 


0.25 


-0.69 


0.03 


51686.632 


21.19 


2.59 


19.51 


2.50 


-1.13 


0.22 


0.59 


0.39 


-0.52 


0.04 


51687.648 


22.97 


1.77 


19.71 


2.94 


-0.88 


0.21 


0.86 


0.29 


-0.39 


0.03 


51688.640 


9.82 


0.31 


11.37 


3.17 


-2.34 


0.35 


-0.06 


0.43 


-0.48 


0.07 


51913.860 


13.86 


2.13 


13.69 


1.88 


-1.46 


0.41 


0.55 


0.19 


-0.23 


0.04 


51918.775 


15.65 


2.15 


15.90 


0.43 


-1.86 


0.29 


0.15 


0.23 


-0.22 


0.04 


51923.824 


14.89 


0.56 


14.91 


0.54 


-1.12 


0.40 


-0.20 


0.17 


-0.14 


0.01 


52019.578 


16.00 


1.81 


16.07 


0.97 


-23.82 


0.45 


-4.06 


1.09 


-1.83 


0.02 


52026.586 


14.48 


1.12 


15.03 


1.05 


-3.80 


0.40 


-0.26 


0.26 


-0.55 


0.05 


52031.573 


14.43 


0.28 


14.64 


0.95 


-9.18 


0.68 


-1.45 


0.30 


-0.64 


0.04 


52319.745 


12.63 


0.36 


12.66 


1.40 


0.26 


0.31 


1.22 


0.12 


-0.26 


0.04 


52322.742''^ 


27.48 


2.56 


19.18 


3.41 


-8.56 


0.40 


-1.23 


0.75 


-1.02 


0.06 


52372.606 


14.70 




9.92 


2.40 


-29.09 


0.95 


-5.60 


0.55 


-2.33 


0.18 


52385.529 


16.44 


0.87 


16.54 


0.62 


-3.82 


0.32 


-0.60 


0.31 


-0.49 


0.02 


52395.617 


4.33 


1.22 


9.46 


2.23 


-20.92 


1.11 


-5.22 


0.71 


-1.97 


0.11 


52396.614 


9.93 




16.20 


3.49 


-11.51 


1.29 


-0.50 


0.50 


-1.78 


0.23 


52710.601 


21.73 


0.30 


20.00 


1.93 


-0.93 


0.22 


0.57 


0.16 


-0.31 


0.03 


52717.620 


14.13 


1.51 


15.79 


1.78 


-13.51 


0.45 


-2.06 


0.50 


-1.16 


0.07 


52723.588 


15.96 


2.33 


15.80 


0.51 


-2.09 


0.31 


0.62 


0.47 


-0.35 


0.01 


53097.541 


15.57 


0.87 


19.17 


1.55 


-6.14 


0.34 


-0.18 


0.25 


-0.72 


0.02 


53104.652 


11.28 


1.18 


11.38 


2.22 


-0.93 


0.21 


0.77 


0.15 


-0.20 


0.10 


53122.685 


15.93 


0.48 


11.70 


3.40 


-22.39 


0.50 


-3.59 


1.23 


-2.00 


0.27 


53137.672 


22.23 


0.66 


15.42 


0.86 


-4.21 


0.21 


0.40 


0.29 


-0.23 


0.11 


53501.652 


29.06 


0.87 


16.43 


1.11 


-8.22 


0.34 


-1.23 


0.67 


-0.59 


0.06 


53524.597 


6.63 


1.35 


7.95 


2.57 


-5.55 


0.27 


-0.18 


0.41 


-0.51 


0.02 


53552.546 


14.31 


1.09 


14.36 


0.37 


-1.06 


0.25 


-0.04 


0.23 


-0.22 


0.04 


53602.488 


29.51 


1.95 


20.32 


4.58 


-8.66 


0.34 


-0.77 


0.60 


-0.68 


0.05 



t He I 5876 A detected 



from the CCF analysis described in Sect. |4.1[ The template 
spectrum was broadened artificially by the v sin i value of 
37kms~^ and then subtracted from each of the TCha 
spectra. We assumed that the photospheric spectrum of a 



TTS resembles that of a main-sequence star of the same 
spectral type. The photospheric subtraction method also 
allows us to verify the possible presence of "veiling". No 
sign of continuum excess emission was found in T Cha 
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Fig. 5. Spectrum of T Cha obtained on 14 Feb 2002. The green line represents the spectrum of the template (HD152391 
G8.5V) artificially broadened for rotation. In each panel is also shown the residual spectrum after subtraction. Slight 
residuals in strong spectral features are mainly due to uncertainties in the normalization to the continuum. The 
spectrum, with a good S/N, shows a redshifted absorption component at H(5. 



(see Fig. 20), confirming the absence of UV-excess in the 
broad-band photometry (see Sect. |8] and Covino et al. 
19961. 



5.3. Emission line spectrum 



Examples of the residual profiles of the main emis- 
sion lines at different dates are shown in Fig. [6] Some 
lines show absorption components shifted toward shorter 
and/or longer wavelengths. The Ha line is the most vari- 
able in both intensity and shape. The residuals of the 
NaiD lines show either redshifted or blue-shifted absorp- 
tion components (that are easily distinguishable from the 
sharp, undisplaced, interstellar component), apparently 
linked to the presence of a correspondent feature in Ha 
and H/3. The sharp interstellar absorption components of 
the D lines from the associated cloud are also distinguish- 
able, when not suppressed by the sky emission. 



Gregorio-Hetem et al. (1992), Alcala et al. (1993) and 



Covino et al. ( 1996 1 already pointed out some outstanding 
properties of the variability in the emission lines, but 
here we exploit the wide spectral coverage of FEROS 
to analyse more lines simultaneously and identify the 
line profile changes at higher resolution than in previous 
works. 



5.3.1. The Ha line 

The most prominent emission line in the spectrum of 
TCha is Ha, which also exhibits impressive variability. 
The observed equivalent width of the line measured in 
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the FEROS spectr£0 ranges from about 0.3 A (on 14 Feb 
2002) to about -30 A (on 8 Apr 2002;^ Besides the vari- 
ations in strength, from pure emission to absorption, the 
line profile is also highly variable in its structure. However, 
we do not observe any correlation between the profile 
of the line and its EW. The line extends approximately 
from —300 to -1-300 km s^^. However, much of the emis- 
sion is concentrated in an interval of almost 400 km s^^ 
in width centred on zero velocity. In most of the spec- 
tra, the Ha line exhibits both blue-shifted and redshifted 
absorption components superimposed on a nearly centred 
symmetric broad emission. The blue absorption is usu- 
ally weaker than the red one, although in three spectra 
of January 2001 it appears stronger (e.g., date 51918 in 
Fig. [6|. During five different epochs, the blue absorption 
reached below the continuum producing a P Cygni pro- 
file. In most cases, when the emission component was not 
so strong, the red absorption was below the continuum, 
producing a typical inverse-P Cygni profile. 

The red absorption is almost always present but with 
different velocity displacements relative to the rest wave- 
length, whereas the blue absorption sometimes becomes 
barely visible or disappears completely. In contrast, the 
emission on the blue side is generally more intense than 
the red one. This suggests that the red wing of the 
line is affected by significant absorption that causes the 
emission to appear weaker than in the blue wing. No 
clear regularity appears to exist in the line changes, but 
some profiles reappear after some time. 

The lack of both detectable veiling and UV excess in- 
dicates that the accretion rate is rather low. Assuming 
that the Ha emission is entirely due to accretion, we es- 
timated the mass acc retion rate. Mamr. by adopting the 



Natta et al. 



(20041 between Mac 



relationship found by 
and the width of the Ha line at 10% intensity. We em- 
phasise that this quantity does not represent an accu- 
rate mass-accretion estimator for two main reasons: i) it 
is not easy to measure, mainly due to the difficulty in 
defining the continuum level, and ii) because the rela- 
tionship, derived for the substellar mass regime, shows 
a larger dispersion at higher masses. We determined the 
width at 10% intensity of the Ha emission only in spectra 
for which we could perform a Gaussian decomposition of 
the line profile into an emission component with overlap- 
ping absorptions. From those, we derived a mean value of 
Maecr=4x lO-^Moyr-i. 



5.3.2. The H/3 line and higher Balmer lines 

The strength of H/3 is also variable and its EW correlates 
well with that of Ha, as shown in Fig. [7j The range of vari- 
ation is between about 1.2 A (on 14 Feb 2002) and about 
—6 A (on 8 Apr 2002). In most cases, the H/3 line presents 
typical double-peak emission with higher emission in the 



blue peak than in the red one. Apart from when emission 
is stronger, the H/3 profile resembles the Ha one. However, 
the inverse P Cygni profile appears more frequently than 
in Ha. 

Even if Ha and H/3 have different optical thicknesses 
and formation regions, their EWs are strongly correlated 
with each other (Fig. [t] lower panel). 

Because of the lower efficiency of FEROS at shorter 
wavelengths and the high variability of TCha in the U 
and B bands, little can be said about the behaviour of 
the higher members of the Balmer series, since H7 and 
US are discernible only in a few spectra. Almost no emis- 
sion is apparent in the residual spectrum, except in one or 
two cases, but the still low S/N prevents us from drawing 
any robust conclusion. The most intense absorptions, ei- 
ther blue- or redshifted, depending on epoch, can also be 
identified in the higher Balmer lines (cf. US in Fig. [5]). 

5.3.3. Forbidden Oxygen lines 

The forbidden [Oi] 6300 and 6363 A emission lines are ob- 
served in the spectrum of T Cha. The contribution of the 
sky emission lines to the FEROS spectra was subtracted 
after scaling the sky spectrum by a factor accounting for 
the different efficiencies of the object and sky fibers, de- 
termined from the ratio of the flat-field intensities in the 
corresponding fibers. Figure |6] shows the sky-subtracted 
line at 6300 A. The width of ^^lOOkms^^ places the line 
just on the edge between low- and high-velocity compo- 
nents for T Tauri stars ( Hartigan et al.|fl995 l. 

As in the case of the Balmer lines, the forbidden oxygen 
lines are of variable intensity, but their profiles remain 
about the same. The most frequently observed profile is 
a symmetrical and well-centred emission, but sometimes 
the line is slightly redshifted by about 10-20 km s^^ (see 
Fig. |6|. The symmetrical profile suggests that both the 
approaching and receding parts of the flow contribute to 
the observed line, i.e., the star-disc system is viewed at 
relatively high inclination. Redshifted [Oi] 6300 A profiles 
have also been observed in BP Tau and RW Aur ( Hartigan 
et al. 19951, but are rather unusual for T Tauri stars. 



Some authors explain the redshifted component in terms 



of asymmetric outfiows (e.g., Hirth et al. |1994| Hartigan 
eraI1[l995l ) 



Remarkably, the equivalent width of the [Oi] 6300 A is 
correlated with that of the Ha emission (see Fig. [t] upper 
panel). As discussed in the following sections, this cor- 
relation indicates that the variations must be related to 
changes in the intensity of the underlying stellar contin- 
uum due to variable circumstellar extinction. 

To verify the hypothesis of variable obscuration, we 
used the only simultaneous spectroscopy and photometry 



data available to us (Covino et al. 19961 to evaluate the 



luminosity in the [Oi] line. 



reported in Tableji] only available in electronic form 
^ conventionally, a negative EW indicates line in emission. 



L[oi] Wioi] ■ /r ■ 47rd 
Wioi] ■ 10-°-^ " ■ 2.69 X lO^^erg s-\ 
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RR Tauri (Rodgers et al. 2002 ). Therefore, the strong vari- 



Fig. 7. Correlations between intensities of the most 
prominent emission lines measured in the FEROS spec- 
tra. 

Table 3. Equival ent widths of [Oi] 6300A measured in 
CASPEC spectra ( Covino et al.|1996 l and line luminosity 
derived using simultaneous R-band photometry. 



JD 


R 


W[0/] 


L[o/] 


-2440000 


mag 


(A) 


(lO^^ergs-i) 


9384.62 


11.47 


-0.78±0.06 


5.4±0.4 


9384.70 


11.70 


-0.91±0.08 


5.1±0.4 


9384.87 


11.88 


-1.24±0.16 


5.9±0.8 


9385.62 


10.56 


-0.32±0.03 


5.1±0.5 


9385.76 


10.41 


-0.25±0.03 


4.6±0.5 


9385.87 


10.39 


-0.24±0.03 


4.5±0.6 



where W[o i] is the measured equivalent width of [Oi] , d 
is the distance of 100 pc, and fu is the observed flux of 
the star in the R-band. The latter is given by fn = ■ 
^0-o.4.fl^ where = 2.26 x lO'^erg s '^cm-^A-M s the 
R-band flux of a zero-magnitude star (Bessel 1979). The 
measured EWs and relative luminosities are reported in 
TableEl 

The luminosity (or the flux) in the [Oi] line remains rel- 
atively constant with time, as in the case of the UXOR star 



ability in the emission-line intensity is caused by changes 
in the stellar continuum flux level rather than the intrin- 
sic variability of the line-emitting region. This explains the 



correlation between Ha and [Oi] discussed in Sect. 5.3.3 
justifying the fact that the intensities of the two lines, 
originating in distinct regions, respond simultaneously to 
the brightness variations of the star. As the analysis of 
flux-calibrated low-resolution spectra presented in Sect.|7] 
shows, the intensities of the Ha and [Oi] lines are also 
found to be related directly to the amount of circumstel- 
lar extinction affecting the photospheric continuum. 

5.3.4. Other emission lines 

Franchini et al. (1992) detected a weak variable emission 
component in the core of the NaiD lines. In our spectral 
subtraction analysis, we found no sign of this emission, but 
instead, were able to identify strong redshifted absorption 
in many residual spectra, that is visible only when the 
redshifted absorption in Ha is also present (although the 
reverse is untrue), independently of the equivalent width 
of the Ha line. On three dates, namely 9 Jan 2001, 21 Apr 
2002, and 9 Apr 2005, when the blue- shifted absorption 
component in Ha was at least as strong as the red one, 
the feature appeared blue-shifted, but there was no sign 
of absorption on 14 Jan 2001, when Ha showed a clear 
P Cygni profile. 

The [Nil] 6583 A forbidden emission line is occasion- 
ally observed in the low-resolution spectra when the star 
reaches its faintest state (see Sect.jT]). 

The He I emission line at 5876 Awas seen on only two 
dates (19 May 1999 and 17 Feb 2002), while missing in 
other spectra close to those dates. This line is probably 
caused by episodic, short-lived flare activity, common in 
young solar-type stars. As shown in Fig.|6] the line appears 
blue-shifted by about — 30kms~^. 

Finally, the emission cores of the two lines 8498 and 
8662 A of the Can IR triplet are clearly visible. The other 
line, at 8541 A, falls in a wavelength gap in-between the 
echelle orders of FEROS. The equivalent widths of the two 
measurable lines remain almost constant, with small vari- 
ations presumably related to chromospheric activity. No 
correlation is found with other emission lines. The pro- 
files appear always symmetrical and well centred on the 
star velocity, indicating presumably a chromospheric ori- 
gin, typical of young stars. 

As for higher Balmer lines, the low S/N prevented us 
from analysing the Can H and K lines, but in some spectra 
of sufficiently high S/N, a strong emission core, indicative 
of intense chromospheric activity, is observed. 

6. Line profile variability on short timescale 

The large variety of Ha profiles shown by T Cha appear 
to occur rather erratically, although this may be caused 
at least partly by the quite uneven temporal sampling of 
our spectra. 
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Figure [9] shows the trends of the equivalent widths of 
the main emission hues with time. Ha and H/3 mainly vary 
in similar ways. On two dates, their behaviour is different, 
i.e., the last day of the first run, when the Ha intensity 
increased as the H/3 intensity decreased, and the third day 
of the second run, when an increase in Ha intensity did 
not correspond to an increase in that of H/3. The [Oi] line 
trend is similar to that of Ha except for during the fourth 
night of the first run and the last night of the third run, 
when both Balmer lines strengthened in emission, while 
the [Oi] lines did not change. The similarities in the vari- 
abihty of these hues indicate that the processes producing 
the observed changes must affect contemporaneously cir- 
cumstellar zones located at different distances from the 
star, i.e., on length scales of between a few tenths of AU 
for Ha and a few AUs for [Oi]. 



Fig. 10. Variations in position of the central peak emis- 
sion and the two absorptions of Ha line during the period 
from 27 to 30 March 1999. The scale on the right refers 
to the red absorption, while the one on left refers to the 
central emission and the blue-shifted absorption. 



We used three different series of spectra acquired with 
a time sampling shorter than one day to investigate vari- 
ations on short timescales. We analysed separately each 
of the following three runs: 1) from 27 March 1999 to 1 
April 1999 (total of 16 spectra in 5 nights); 2) from 16 
to 22 May 1999 (total of 6 spectra in 6 nights); and 3) 
from 9 to 23 May 2000 (total of 5 spectra in 5 nights). 
Figure [s] shows in each panel the overlap of Ha (top-left 
box) and H/3 (bottom-left box) line profiles for each of 
the three runs. The average normalized profile (/„(A)) is 
also superposed, and is represented by a thicker line. The 
profiles, and their behaviour, appear to differ in each run, 
being more erratic in the second period than in the other 
two. The Ha and H/3 line intensities and profiles change 
on a daily timescale, with only minor changes occurring 
during the night. We report in Fig. |10[ the variability in 
the positions of the central peak and the two absorption 
components in Ha during the first run. 

To quantify the variability in each line, we examined 
the variance, cr, and the normalized variance, tT„, in the 
different runs, defined respectively as: 



Sf.i[/„,.(A)-(/„(A))]7(iV~l) 



(1) 



and 



aUX)^<j\X)/{UX)) (2) 

where /„^i(A) is the i-th spectrum. These quantities are 
plotted in the right panels of Fig.|8] If variability were 
merely caused by the brightening or fading of the stellar 
continuum flux, the shape of normalized variance profiles 



would be the same as the average line profile (Johns & 



6.1. Correlation Matrices 

To test whether the variations across the emission line 
profiles have a common origin, we computed their correla- 
tion matrices (CMs). CMs indicate the linear correlation 
coefficient between the variation in each velocity bin of 
the spectral line profile and variations in all other bins for 
the same or for two different lines. A strong correlation 
between different velocity bins is indicative of a common 
origin, or emitting region. If the variability in the emission 
line is linked to variability in the continuum, then the line 
profile is highly correlated over a wide range of velocity, 



Basri 1995 1. This is generally not the case, suggesting that 



some other changes in the Ha forming region also occur. 



producing a typical squarish shape in the CM (Johns & 
|Basri„ 19951 |Alencar fc Batalha„2002) . 

We computed the CMs of Ha and H/3 with respect to 
themselves and to each other. Two-dimensional plots of 
these matrices for Ha and H/3 in three different runs are 
shown in Fig. |11[ No anticorrelation is evident in the ma- 
trices. In each run, there are regions that correlate well 
with themselves, but do not correlate with the remainder 
of the line. No significant correlation is seen between the 
red (> 200kms~^) and the blue (< — 200kms~^) wings 
of the profile. However, in the second run, the Ha line ex- 
hibits a correlation over a wide velocity range, its squarish 
form extending from about —300 to nearly 150kms~^. 
The correlation is disrupted by the blue-shifted absorp- 
tion, from —150 to — 50kms^^, and the redshifted com- 
ponent, up to about — lOOkms"^, which do not corre- 
late with the remainder of the line profile but with the 
redshifted emission peak, around 130kms^^. This be- 
haviour is not observed in the other epochs. Similarly, 
in the first run, the broad red absorption from —20 to 
150 km s^^ breaks the correlation with the blue emission 
wing, which still correlates weakly with the red emission 
wing. Interestingly, the square regions in the CMs corre- 
spond to some of the highest peaks in the variance profiles, 
i.e., the blue emission wing, the central peak and the red 
absorption, respectively. 

The H/3 CMs are more affected by noise due to lower 
S/N at those wavelengths, but again traces of a squarish 
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form can be seen in particular in data from the second 
run. We also computed the CM between Ha and H/?, but, 
although the profiles of the two lines appear to be sim- 
ilar, the corresponding correlation matrices are affected 
too significantly by noise and no useful information can 
be extracted from them. 

In conclusion, coherent changes in the velocity inter- 
val between —300 and 120 km s^^ appear consistent with 
variations dominated by variable circumstellar obscura- 
tion of the stellar photosphere, giving rise, due to a con- 
trast effect, to changes in the relative intensity of emission 
lines. In the case of magnetospheric accretion, parts of the 
Ha profile, i.e., the redshifted component, are expected to 
originate in regions close to the star and eventually to 
be affected by direct obscuration. These may be identi- 
fied kinematically with the intervals where the correlation 
breaks down whereas the differences from one epoch to 
the other possibly reflect changes in the geometry of the 
emitting region. 

7. Variability in the low-resolution spectra 

Strong variability was also detected in our low-resolution 
spectroscopy. We refer in particular to Table where we 
report the equivalent widths of the Ha and [Oi]A6300A 
lines, and to Fig. [12] where a sample of low-resolution 
spectra is showrj^ 

We note that during the May-June 1993 run, the star 
exhibited spectacular changes between one night and the 
next. For instance. Ha was observed to vary between weak 
absorption (EW« 0) on May 31, to a very strong emission 
line (EW~ — 50A) on June 1. The latter spectrum was also 
characterised by strong [Oi] lines, and the [Nil] line was 
also present, while no emission lines were detected in data 
acquired the two nights before. 

The spectacular variability of these lines is illustrated 



Table 4. Equivalent widths of Ha and [Oi] 6300A mea- 
sured on the low-resolution spectra (see the atlas in elec- 
tronic form) and corresponding total extinction. Ay, ob- 
tained as described in Sect. [8] Only in electronic form. 



in Fig. 13 where two low-resolution spectra, correspond- 
ing to close to the maximum and minimum brightness of 
the star, just two nights apart, are compared after nor- 
malising each of them to the flux at A = 5880 A. We re- 
call that these spectra are calibrated in relative flux (see 
Sect. 3.2 1. The variation in the continuum slope between 
the two dates is also remarkable, in addition to the emis- 
sion features appearing prominent when the photospheric 
continuum looks fainter and heavily reddened. 

However, no change in the spectral type of the star is 
observed. To verify this, we used low-resolution spectra 
calibrated in relative flux, and applied the methods de- 



scribed in Alcala et al. (20061 and Gandolfi et al. (20081 



to determine simultaneously the spectral type and the 
visual extinction at the phases of minimum and maxi- 
mum brightness, respectively. The result, shown in Fig. |14| 
is that the spectral type did not vary to within about 
half a sub-class, while the visual extinction changed from 

® available only in electronic form 

^ the full atlas of low-resolution spectra is reported in 
Fig. fig:lowres-atlas, only in electronic form 



Date 



JD 
-2440000 



Wj 



(A) 



Oh a 

(A) 



w 



[Oi_ 



(A) (A) 



07 May 92 

08 May 92 

10 May 92 

11 May 92 

12 May 92 
27 Mar 93 

27 Mar 93 

28 Mar 93 

28 Mar 93 

29 Mar 93 

29 Mar 93 

30 Mar 93 

30 Mar 93 

31 Mar 93 
01 Apr 93 

01 Apr 93 

02 Apr 93 

02 Apr 93 

30 May 93 

31 May 93 
31 May 93 

01 Jun 93 

03 Jun 93 

04 Jun 93 

05 Jun 93 

05 Jun 93 

06 Jun 93 

28 Jun 93 

29 Jun 93 

02 Jul 93 

03 Jul 93 
24 Mar 94 

24 Mar 94 

25 Mar 94 

25 Mar 94 

26 Mar 94 

26 Mar 94 

27 Mar 94 

27 Mar 94 

28 Mar 94 

28 Mar 94 

29 Mar 94 
29 Mar 94 

09 Apr 95 

10 Apr 95 

11 Apr 95 

12 Apr 95 



8749.668 
8750.606 
8752.628 
8753.634 
8754.600 
9073.698 
9073.847 
9074.689 
9074.843 
9075.680 
9075.873 
9076.758 
9076.840 
9077.693 
9078.700 
9078.879 
9079.687 
9079.879 
9137.531 
9138.515 
9138.738 
9139.545 
9142.495 
9143.494 
9143.736 
9144.496 
9144.718 
9167.463 
9168.478 
9170.519 
9172.467 
9435.745 
9435.887 
9436.698 
9436.900 
9437.532 
9437.672 
9438.519 
9438.662 
9439.520 
9439.661 
9440.518 
9440.661 
9816.807 
9817.789 
9818.805 
9819.788 



-4.58 
-3.63 
0.40 
-4.87 
-6.42 
-1.09 
-0.51 
1.80 
2.08 
0.28 
-0.14 
-23.70 
-34.00 
-9.20 
-1.80 
-0.32 
1.90 
1.64 
1.40 
0.26 
-0.25 
-44.60 
0.16 
-0.23 
0.66 
-6.20 
-5.81 
-2.40 
-19.3 
-7.21 
-57.0 
1.73 
2.20 
0.55 
-1.03 
0.96 
0.21 
-1.10 
-0.95 
1.61 
0.95 
0.25 
1.05 
-0.19 
-38.70 
0.70 
1.05 



0.30 
0.40 
0.10 
0.70 
0.80 
0.20 
0.20 
0.20 
0.20 
0.20 
0.10 
2.00 
4.00 
1.00 
0.10 
0.10 
0.10 
0.10 
0.10 
0.10 
0.10 
5.00 
0.10 
0.10 
0.10 
0.50 
0.60 
0.30 
2.00 
0.60 
5.50 
0.10 
0.20 
0.10 
0.10 
0.10 
0.10 
0.10 
0.10 
0.20 
0.10 
0.10 
0.20 
0.10 
3.50 
0.10 
0.10 



-0.33 
-0.35 
-0.13 
-0.22 
-0.87 
-0.23 
-0.19 

-0.06 
-0.25 
-0.11 
-4.16 
-6.29 
-1.27 
-0.26 
-0.14 
-0.06 



-0.26 
-0.35 
-9.20 
-0.36 
-0.48 
-0.23 
-0.57 
-0.68 
-0.35 
-4.18 
-0.93 
-9.00 
-0.09 
-0.10 
-0.23 
-0.36 
-0.24 
-0.36 
-0.30 
-0.27 



-0.03 
-0.03 
-0.04 
-4.31 
-0.05 
-0.14 



0.05 
0.05 
0.03 
0.05 
0.10 
0.05 
0.03 

0.05 
0.05 
0.03 
0.50 
0.55 
0.15 
0.05 
0.05 
0.05 



0.05 
0.05 
0.55 
0.05 
0.05 
0.05 
0.05 
0.05 
0.05 
0.50 
0.10 
0.60 
0.05 
0.03 
0.05 
0.05 
0.05 
0.05 
0.05 
0.12 



0.05 
0.05 
0.05 
0.45 
0.05 
0.05 



Av 
(mag) 



2.2 
2.4 
1.8 
2.0 
2.8 
1.7 
1.4 
1.5 
1.5 
1.7 
1.7 
3.3 
3.6 
2.6 
1.9 
1.3 
1.5 
1.2 
1.2 
2.0 
2.4 
4.6 
1.5 
2.3 
2.2 
2.8 
3.2 
1.8 
3.1 
2.1 

X 

1.8 
2.1 
1.5 
1.6 
1.7 
1.5 
1.8 
1.6 
1.3 
1.7 
1.2 
1.6 
1.6 
3.4 
1.5 
1.8 



X Ay not determined due to very faint continuum 



1.2 mag in the bright phase to about 4.6 mag in the faint 
state, respectively. This corresponds to an extinction in- 
crease of at least 3.4 mag between the bright and faint 
states. By instead using the spectrum corresponding to 
the bright level as a template, an extinction of 3.3 mag is 
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5000 



Fig. 13. Low-resolution spectra of T Cha near minimum 
(June 1, in red) and maximum brightness (May 30, in 
blue). The spectra are caHbrated in relative flux and nor- 
malised to the flux at A = 5880 A. Emission features, only 
visible in the spectrum corresponding to the fainter level, 
are marked. 



needed to reproduce the faint spectrum, and, apart from ^ 



emission lines, the residual is indeed quite low (cf. Fig. 15 ) 



We applied the same procedure to each of the low- 
resolution spectra and determined the corresponding value 
of Ay (see Table |4]), with estimated errors of about 10%. 
This takes account of the uncertainty in the fit, as well 
as the relative flux calibration, but not systematic effects 
caused by deviation of the circumstellar extinction from 
the normal interstellar law, as shown in Sect. [8) The two 
panels of Fig. [16] show that a clear trend exists between 
the amount of visual extinction and the intensity of the 
Ha and [Oi] 6300A emission lines. This indicates that the 
observed variations do not reflect intrinsic changes in the 
stellar photosphere, but arise presumably from variable 
circumstellar extinction, as can be inferred in Fig. |13| from 
the different continuum slopes of two spectra taken near 
maximum (on 30 May 1993) and minimum (on 1 Jun 1993) 
brightness. The histogram in Fig. [17] indicates that the 
highest extinction events are relatively rare, while events 
with differential extinction below 1 mag are more frequent. 



8. The spectral energy distribution 

We constructed the observed SED of T Cha using all op- 
tical and near-IR photometry available to us (Alcala et 



al. 1993 Covino et al. 19961, as well as data from public 



catalogues. In Fig. [18j the SED is shown for wavelengths 
shorter than 10 /im. We remark that photometric errors 
are smaller than the symbol size, and the considerable 
scatter of the optical {UBVRI bands) data points reflects 
the significant variability of the star. Although of smaller 
amplitude, variability in the near-IR [JHK bands) can 
also be appreciated. Therefore, relying on data that have 
not been acquired simultaneously may be a limitation in 
the SED analysis. 



1993 June 1 




5500 



6000 
A (A) 



6500 



Fig. 14. Best-fit of spectral type templates to the spectra 



shown in Figure 13 The spectra of T Cha are shown as 
thin (black) lines, while the templates arc shown as thick 
(red) lines. The corresponding values of total visual ex- 
tinction, as derived from the two-parameter fit adopting a 
normal interstellar extinction law with Ry = 3.1, are also 
indicated for both the faint and bright stage, respectively. 



At the brightest level, the observed J flux was lower 
than the expected stellar photospheric flux. However, we 
note that T Cha has not been monitored as extensively in 
the near-IR as in the optical. Therefore, the available data 
probably did not detect the brightest J flux. For that rea- 
son, the J flux was excluded from our SED-fitting. Excess 
emission is instead observed for wavelengths longer than 
the H band. The agreement of the Spitzer spectroscopy 
with the photometric data in the mid- and far-IR (cf. 



Fig. 19 1 indicates that there is no significant variability at 



longer wavelengths. This is expected in the case that the 
strong brightness and colour variations affecting the star 
are due to variable extinction of the stellar photosphere 
from inhomogeneous, circumstellar material. 

By assuming variable extinction, we used relationships 
derived from U BVRI photometry (Covino et al. 1996) to 
probe the dust column affecting the brightness of T Cha. 
In Table [5] we report the circumstellar extinction law for 
T Cha, expressed by the total to selective extinction ratios, 
R\ — A\/ Eb-v , derived from the differential brightness 
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0.5 



0.001 +0.03 



[01] 



[01] 



Ha 



5500 



6000 



6500 



MA) 



Fig. 15. The top panel shows the spectrum at the faint 
state as a thin (black) line, and the spectrum at the bright 
stage, used as a template, is shown as a thick (red) line. 
The lower panel shows the residual. The average residual 
and rms are indicated. The emission lines, not included in 
the deriving the average value of the residual, are marked. 




Fig. 17. Histogram of differential visual extinction. The 
most frequent value is around 0.5 mag, while the maximum 
observed is 3.3 mag. 



Table 5. Total to selective extinction ratios, R\ — 
A\/Eb-v, for TCha, derived from UBVRI photometry 
(Covino et al. 1996). 



Ru 


Rb 


Rv 


Rr 


Ri 


8.4 ±0.4 


6.5 ±0.3 


5.5 ±0.2 


4.9 ±0.2 


4.5 ±0.2 




Fig. 18. Observed SED of TCha for A < l^m. The 



U,B,V,R,I data (filled circles) are from Alcala et al. 
( |1993[ ); [Covino et al] (|1996[). T h e IR fiuxes come from 
difi^erent papers: Alcala et al.| ( 1993[ ) (open squares), 
2MASS fiuxes (open circles), DENIS survey (filled tri- 
angles) , [Morel fc Magnenat Photometric Catalog ( 1978[ ) 
(filled squares), Spitzer IRAC bands (filled diamonds), 
MSX6C (crosses). IRS spectra are also superposed. The 
upper and lower dashed lines represent the SEDs of a 
G8V standard, reddened to the corresponding {B — V) 
values that match the brightest and faintest /-band ob- 
served fiuxes, respectively. 



variations measured in the U BVRI bands (Covino et al. 
1996). The resulting value of Rv deviates from the nor- 
mal value {Rv = 3.1) for the diffuse interstellar medium, 
which indicates different physical properties of the circum- 
stellar grains. In particular, a higher value of Rv means 
a fiatter extinction curve (i.e., grayer extinction), and re- 
veals the presence of larger grains, thus implying probable 
grain growth and depletion of small grains ( |Draine|[2008t . 



Therefore, in the following SED analysis, we consid- 
ered more suitable the fiux measurements corresponding 



to the brightest level of the star, shown in Fig. 19 



Disc models such as those by RobitaUle et al. (2006 



2007) do not take account of the presence of a partially 
evacuated inner hole or gap, so may be inadequate in re- 
producing the SED of transitional objects such as T Cha. 



Brown et al. (20071 found that the mid-IR spectrum of 
TCha can be reproduced by a disc truncated at 0.08 AU 
with a gap between 0.2 and 15 AU, but they did not con- 
sider the strong variability of the object nor use the avail- 
able millimetre data. 
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estimates of the disc mass based on single measurements 
at millimetr e wavelengths are in the range of 3xl0 ~'^- 
IxIO^^Mq dHenning et al. 



1993 



Lommen et al. 



20071 



but our estimate reproduces well the observed data for 
wavelengths longer than 60 /im. The disc radius is more 
poorly constrained, but values between 100 and 150 AU 
provide quite reasonable results. The latter is consistent 



with the 3.3 mm observations by Lommen et al. (20071 



that did not resolve the disc around T Cha. We note that 



the disc radius adopted by Brown et al. ( 2007 1 is a factor 



1 2 
log X(/im) 



Fig. 19. Observed data-points used for SED fitting. 
Optical: symbols as in Fig. [18] IRAS: open triangles. Also 
superposed are IRS spectra (thin line). Three disc models, 
whose parameters are reported in Table [6] are represented 
by the solid, dashed and dotted lines, for different disc 
radii. 

Table 6. In the upper table, we report the estimated stel- 
lar properties, and in the lower table, the stellar and disc 



of two higher than the value derived here. 

Although the model with a disc radius of 150 AU re- 
produces well the observed SED at optical, far-IR, and 
millimeter wavelengths, all the models predict too high 
flux between 10 and 30 /im relative to the observed SED. 
Thus, despite our attempt to minimize the eS'ects of colour 
changes by adopting the brightest flux level in each band, 
we were unable to fit the complete SED of T Cha with a 
simple reprocessing-disc model. This failure is ascribed to 
excess emission in the near-IR. Therefore, the most plau- 
sible disc model for reproducing the SED of TCha is a 



disc with a gap, as proposed by Brown et al. (2007) 



parameters resulting from SED-fitting using , DuUemond et| 9. Discussion and conclusions 



al. (20011 model 



Stellar Parameter 



value 



L (Lq) 

Teff (K) 

R (Ro) 
M (Mq) 
Age (Myr) 



2.& 

5520* 

1.8* 

1.3" 

10" 



Model Disc Best-Fit Parameters 



U (Lq) 

T* (K) 

M. (Mq) 
Mdisa (Mq) 

Rout (AU) 

Rrn (AU) 

Temperature of rims 
Inclination i 



2.5 

5400 

1.5 

5x10"^ 

150 

0.13 

1300 
75° 



f from SED integration and adopting d=100pc and Av = 1.3; 
X from the relations by |Kenyon &: Hartmann| ( |1995[ ) for a G8 V 
star; * from the Stefan-Boltzman law; ** from comparison with 



Palla & Stabler ( 1999 1 PMS tracks 



We explored the possibility of modelling the SED, si- 
multaneously from optical to millimetre wavelengths, us- 



Before discussing our interpretation of T Cha, we briefly 
summarise the most signiflcant properties of the star de- 
rived from analysis of the spectroscopic observations pre- 
sented here, complemented with published photometry 
and other data from public archives: 

a) the spectral type G8 and a lithium abundance indica- 
tive of a very young age are confirmed 

b) RV variations of almost lOkms^^ on relatively short 
timescale are detected, although without revealing any pe- 
riodicity that could be releted to orbital motion. Besides 
that, the v&ini determinations obtained from FEROS ob- 
servations are consistent with the value of 37±2kms~^, 
except for a few spectra in which the photospheric lines are 
significantly narrower, yielding a vsmi of 30±2kms~^. 

c) strong variability in the main emission features (e.g.. 
Ha, H/3, and [Oi] lines) on timescales comparable to that 
of photometric variability, and a close correlation between 
the Ha and [Oi] line intensities are found. 

d) variations in the photospheric continuum slope arc ob- 
served in connection with intensity changes in the main 
emission features (i.e., the line intensities strengthen when 
the star reddens and fades). 

e) clearly evident relationships between brightness and 



ing the CGPLUS prescription by Dullemond et al. (2001 ). 
Starting from the stellar parameters in Table [6] we con- 
structed a grid of SEDs. The results of the SED fitting are 
reported in Table[6]and the best-fit models are overplotted 
on the observed SED in Fig. [19] The dust opacities cal- 



colours based on optical photometry (Alcala et al. 1993 



Covino et al. 1996| have allowed us to derive a circum- 



culated by Laor & Draine (1993), modified to match the 



Beckwith et al. (1990) opacity law at wavelengths longer 
than 100 /xm were used. We note that a blackbody with 
Teff =5400 K better reproduces the optical data. Previous 



stellar extinction law characterised by Ry = 5.5 and infer 
the presence of large circumstellar dust grains, 
f) the observed SED, showing a dip at mid-IR wave- 
lengths, is typical of a transitional disc. SED modelling 
has confirmed the presence of a disc with a gap and a 
high inclination angle, as first indicated by [Brown et al.| 
(|2§07l). 
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Several of the properties outlined above provide clues 
(e.g., items c, d, and e) that a mechanism of variable cir- 
cumstellar extinction is most probably operating in T Cha, 



as in UXOri-like stars (Natta et al. 2000 Rodgers et al. 15 and 20) suggesting very low or no mass accretion 



significantly decrease with respect to the earlier CTTS 
phase. In T Cha, no veiling is apparently detected, not 
even in the stages of strong Ha emission (see Figs, [sj 



2002 1 . The simultaneous brightening of Ha emission and 



oxygen forbidden lines as reddening increases, as well as 
the existence of a tight correlation between brightness and 
colour variations, can be most plausibly explained in terms 
of variable extinction produced by obscuring clumpy ma- 
terial moving across the line of sight to the star. First of 
all, whereas the Hydrogen Balmer lines originate in regions 
closer to the star, the oxygen forbidden lines form instead 
farther away from the star and are hence not expected to 
be sensitive to variations occurring in the inner parts of 
the disc. Moreover, the dominant timescales and ranges 
of the observed variability phenomena indicate that the 
occulting clumps must be small compared to the size of 
the disc and more likely concentrated within a distance 
of a few tenths AU from the star. These clumps occasion- 
ally obscure the star favoured by the high inclination of 
the disc, while not affecting either the outer, extended, 
low-density wind region traced by the [Oi], or the bulk 
of the circumstellar hydrogen emitting zone. The conse- 
quence is a strong contrast enhancement of all emission 
lines relative to the photospheric continuum when the star 
is highly attenuated, while the lines weaken or even dis- 
appear as the photospheric contribution turns brighter. 
Furthermore, these dusty clumps might also be associ- 
ated with patchy gas structures orbiting close to the star 
(< 0.2 AU) that cause daily changes in the Ha line profile, 
similar to those we observe (see Sect. 5.3.11, due to vari- 
ations in the absorbing column of material shielding the 
star. If the Ha variability were merely a consequence of the 
brightening and fading of the stellar continuum flux, the 
shape of the line profile should not change with time. This 
is generally not the case, suggesting that some other mod- 
ifications in the Ha forming region must also intervene. 
In particular, we expect that the parts of the Ha profile 
originating closer to the star (e.g., the redshifted absorp- 
tion components, in the case of magnetospheric accretion) 
also show the effects of occultation by orbiting material as 
well as eventual modulation induced by rotation. This is 
consistent with the nightly line profile changes seen in the 
FEROS high-resolution spectra and the continuum flux 
(i.e. extinction) variations sharing a common timescale. 
However, while we can estimate that the time required 
for an occulting clump to appear in front of the star likely 
ranges between about one-third and two-thirds of a day, it 
is difficult, from the present data, to quantify in a reliable 
way the fraction of Ha emission that may originate close 
to the star, as we are unable to determine information 
about the detailed geometry of the occultation. 

But what is the source for the bulk Ha emission in 
T Cha? It is generally accepted that the strong Ha emis- 
sion line in T Tauri stars is due to mass accretion, which 
leads to UV excess emission and photospheric continuum 



However, the non-detection of veiling in T Cha may also 
caused by the relatively high continuum flux emitted by its 
G8 photosphere. Hence, we cannot exclude the possibility 
that erratic accretion episodes are caused by instabilities 
in the inner disc, in which inf ailing gas produces the vari- 
able redshifted absorption features that we observe in the 
Ha line profile. However, the [Oi] A6300 A emission, which 
corresponds to a luminosity of log(L[oi]/-^0) ~ —4.8, im- 
plies the presence of strong winds, although the width of 
the line (wlOOkms"^, see Fig. [g]) is at the lower limit of 
a typical high- velocity component ( |Hartigan et al.||"l995 l. 
It is thus probable that the bulk of the Ha emission ob- 
served in TCha is also produced in a wind, as reported 
for some Herbig Ae stars (Kraus et al. 2008). Moreover, 



the Ha versus [Oi] correlation shown in Fig. [7 is more eas 
ily explained in terms of variable circumstellar extinction, 
rather than variable accretion. 

By both assuming variable circumstellar extinction, 
and adopting the relation Ay = l.Q^&Cext{V)NH and 
the extinction cross-section per H nucleon for Ry — 5.5, 
C^xtiV) = 6.715 X IQ-'^^cia^/H, given by |Draine| (pOOSt, 



we estimate the column density for an obscuring clump to 
be Nh ~ 1.4 X lO^^Av cm^^. By also assuming the maxi- 
mum observed differential extinction of 3.3 mag, we obtain 
a hydrogen column density Nh = 4.6 x 10^^ cm~^. We can 
estimate the mass of the obscuring clumps, Md, by assum- 
ing that they cover completely the photospheric disc of the 
star. Using the stellar radius reported in Table [6] the re- 
sult is A'hi « 4 X 10^° g, i.e., of the order of 2 x 10"" Mq. 
Similarly, assuming the most frequent differential circum- 
stellar extinction from the histogram in Fig. [17] the result 
would be Mri ~ 3 x 10^^'' Mfr,. In contrast, the mass of as- 



teroid 253Mathildeis5xlO~i^MQ ( [Yeomans et al.|l998 l. 

Finally, could the RV variations in T Cha becaused 
by a low-mass companion? Some PMS stars with transi- 
tional discs, such as CoKuTau/4, have been discovered to 
be binaries (e.g., Ireland fc Kraus| 2008). The inner edge 
of circumbinary discs is believed to be set by dynamical 
truncation. In the case of CoKuTau/4, where the binary 
separation is about 8 AU, the truncation radius would be 
of the order of 13-16 AU, depending on the eccentricity of 
the binary ( Ireland fc Kraus|2008 1 . As concluded from our 
analysis of the CCF bisector, the RV variation of T Cha is 
about lOkms^^. A very low-mass object, say a 0.05 M© 
brown dwarf on a circular orbit with a radius of about 
0.1 AU may cause this variation, producing a disc trun- 
cation radius of about 0.2 AU, consistent with the inner 
disc edge proposed bylBrown et al.l pOOTl. However, no 



veiling (Hartigan et al. 1995). Transitional objects, how- 



ever, have more evolved discs in which accretion should 



obvious periodicity is detected in our RV measurements. 
On the other hand, the gap proposed by |Brown et al.| 
(2007) cannot be explained in terms of truncation by a 
companion at a separation of about 10-15 AU. At such 
separation, a 0.9 Mq star would be necessary to produce 
the observed RV variation. An object of that mass and 
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separation would have been eventually detected either by 



us as a spectroscopic binary, or by Kohler (20011 in his 
speckle survey. Thus, the gap of the transitional disc in 
T Cha is probably not the result of dynamical truncation 
due to a binary companion. In any case, the presence of 
a lower-mass (stellar or substellar) companion cannot be 
definitely excluded because a quasi-periodicity of about 
3 days is clearly evident throughout the entire data-set 
of T Cha. On the other hand, T Cha is a relatively fast 
rotator (vsini 37kms^^) and a magnetically active 
{logLx ~ 30ergs~^) PMS star. Therefore, the observed 
RV variations might be alternatively induced by stellar 
spots, as in the case of the WTTs V410Tau ( [Stelzer et 
al. 2003). Otherwise, we cannot exclude the possibility 



that the same circumstellar inhomogeneities responsible 
for the variable extinction may also alter the symmetry of 
the photospheric lines, and cause apparent RV shifts due 
to a sort of Rossiter effect, as in eclipsing binaries, which 
occurs when the stellar photosphere is partially obscured. 
In all cases, the presence of an inner disc characterised by 
a clumpy structure and the evidence of grain growth both 
indicate that T Cha is a likely candidate to host planet 
formation. 
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Fig. 6. Examples of residual line profiles of T Cha in velocity scale. The vertical dotted line marks the rest position 
in the star frame. 
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Fig. 8. Line profiles {left panels) and corresponding vari- 
ance and normalized variance profiles (right panels) for Ha 
and H/3 in three different runs. 



Fig. 9. Equivalent widths of the main emission lines 
during the runs displayed in Fig. [8] The left scale 
refers to Ha, while the right one is for H/3 and [Oi] 
lines. 



Schisano et al.: Variability of TCha 



21 




Fig. 11. Two-dimensional maps of the correlation matrices of Ha (left panels) and H/3 (right panels) line profiles with 
themselves, for the three runs of Figs. |8] and [9] Different squarish regions are distinguishable in the first two runs, 
indicating a coherence in the variations of Ha intensity in the velocity interval from —300 to 120 km s^^. 




Fig. 12. Low-resolution spectra of TCha. The panels on the left side display the whole range, those on the right 
show the interval containing Ha and the Lithium line. The complete atlas of low-resolution spectra is only available 
electronically. 



5 - 



4 - 



q 1 1 — I I I I I 1 1 



"T 1 — I I I I I 1 1 



cm 

(0 o 



> 
< 



I I I I I L_ 



TTT| 1 1 — I I I I I I I 1 1 — I I I I I I I 1 1 — I I I I I I I 




1 10 
WH„(TCha) - WHa(G8) [A] 



I 1 1 1 1 i Mill I I I I I I I I I 

100 0.01 0.1 1 10 

W[oi] [A] 



Fig. 16. Visual extinction derived from low-resolution flux-calibrated spectra versus the equivalent width of Ha 
(left panel) and [Oi] 6300A (right panel) emission. For Ha, the abscissa refers to the net emission after subtracting a 
constant value corresponding to the photospheric absorption for a G 8 star. The horizontal dashed line indicates the 
minimum value of Ay obtained for T Cha at maximum brightness. 
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Fig. 20. Comparison between T Cha spectra observed in two different epochs: in blue the spectrum observed on f 4 
Feb 2002, in red the one observed on 26 Apr 2004, in black the difference between the two. Various spectral lines at 
different wavelengths are marked. 
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